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Abstract. 

We have obtained empirical relations between the p-mode frequency shift and the 
change in solar activity indices. The empirical relations are determined on the basis 
of frequencies obtained from BBSO and GONG stations during solar cycle 22. These 
£ — ■ relations are applied to estimate the change in mean frequency for the cycle 21 and 

f\j ' 23. A remarkable agreement between the calculated and observed frequency shifts for 

the ascending phase of cycle 23, indicates that the derived relations are independent 
of epoch and do not change significantly from cycle to cycle. We propose that these 
relations could be used to estimate the shift in p-mode frequencies for past, present 
and future solar activity cycles, if the solar activity index is known. The maximum 
frequency shift for cycle 23 is estimated to be 265 ± 90 nHz, corresponding to the 
predicted maximum smoothed sunspot number 118.1 ± 35. 
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■ 1. Introduction 

It is now well established that the solar p-mode oscillation frequencies 
vary with solar activity. The first evidence of this effect came from 
the analysis of low degree acoustic frequencies derived from solar ir- 
radiance data by the Active Cavity Radiometer instrument on board 
Solar Maximum Mission satellite (Woodard and Noyes, 1985). Using 
the Doppler velocity data, Palle, Regulo, and Roca Cortes (1989) also 
showed that the frequency shifts are well correlated with solar activity 
cycle and obtained a shift of 0.44 ± 0.06 fjHz between the minimum and 
maximum of solar cycle 21. Later, several other authors have extended 
these studies to different epochs with new and improved data sets for 
intermediate (Woodard et al., 1991; Bachmann and Brown, 1993) and 
low degree modes (Jimenez- Reyes et al., 1998 and references therein). 
A consistent and continuous data set of intermediate degree p-mode 
frequencies is recently made available from the Global Oscillation Net- 
work Group (GONG) for the period May 1995 to October 1998, with 
an interval of 108 days. Using some of these data sets between August 
1995 to August 1997, Bhatnagar, Jain, and Tripathy (1999) obtained a 
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decrease of 0.06 fiHz in the mean frequency during the descending phase 
of the solar cycle 22 and an increase of 0.04 /iHz during the ascending 
phase of the cycle 23. Howe, Komm, and Hill (1999), using a subset 
of GONG frequencies, also demonstrated that the p-mode frequencies 
vary with solar activity cycle. 

The prime motivation of this work is to derive empirical relations 
between the shift in frequency and change in the level of the activity 
indices. The derived relations could be used to estimate the frequency 
shift for past and future solar cycles. In this paper, we have used the 
change in activity indices corresponding to the same epoch as the fre- 
quency shifts instead of the actual value of activity. We find that the 
correlation between change in activity and frequency shifts is better as 
compared when the actual value of the activity indices were used. Using 
a similar approach Rhodes et al. (1993) showed that the frequency 
shifts between 1981 and 1989 are correlated with the change in various 
activity indices; e.g., the sunspot number, sunspot area, and irradiance 
measurements. 



2. Observational data 



The mode frequencies for cycle 22 in the intermediate degree range 
are available from the Big Bear Solar observatory (BBSO), LOWL 
instrument, South Pole expeditions and GONG project. As shown by 
Jain et al. (1999), frequencies derived from South Pole observations are 
systematically higher as compared with other three data sets. Similarly 
frequencies derived from LOWL instrument are from one year power 
spectra, which may not be valid for the study of short period solar cycle 
variation. Thus, in this study we use p-mode frequencies obtained from 
BBSO and GONG, for the period 1986 to 1996; these are summarised 
in Table I. It may be noted that we have used only the continuous and 
independent data sets from the GONG network. As it was shown by 
Libbrecht and Woodard (1990) that the mode frequency shift strongly 
depends on the frequencies, we have used only the common modes 
available between GONG and BBSO data sets. This selection criteria 
generated a total number of 412 common modes in the frequency range 
of 1500 and 3500 [iHz and i between 5 to 99. Further, we have cate- 
gorised the data sets into two groups according to the data source as 
shown in Table I. 
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3. Analysis and results 

The mean frequency shift is calculated by taking simple difference be- 
tween any two data sets chosen from the same group. With 9 frequency 
data sets (see Table I), the possible combinations (choosing any two at 
a time, out of the total data sets from an individual group) produced 
16 values of the mean frequency shifts. In order to investigate how the 
frequency shifts are correlated with the change in the level of activity 
indices and to derive empirical relations between them, we have used 
seven different activity indices representing the magnetic and radiative 
indices. These indices are: R/, unsmoothed International sunspot num- 
ber obtained from the Solar Geophysical Data (SGD); KPMI, Kitt Peak 
Magnetic Index from Kitt peak full disk magnetograms; MPSI, Mag- 
netic Plage Strength Index from Mount Wilson magnetograms (Ulrich , 
1991); FI, total flare index from SGD and Atag (1999); He I, equivalent 
width of Helium I 10830A line (Harvey 1984), averaged over the whole 
disk from Kitt peak; Fio, integrated radio flux at 10.7 cm from SGD, 
and R s , smoothed International sunspot number obtained from SGD. 
A mean value for each activity index was computed, for the epoch 
corresponding to the actual frequency interval. 

To study the relative variation in the mean frequency shift bv with 
the change in activity index 5i, we assume a linear relationship of the 
form: 

5v = a 5i + b, (1) 

where slope a and intercept b are obtained by performing a linear least 
square fit and is shown in Figure 1. The solid line represents the best 
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Figure 1. Linear regression of mean frequency shift with change in activity indices. 
The solid line is the least square fit. The triangles represent BBSO data while squares 
are for GONG data. 



regression fit and confirms that the data sets are consistent with the 
assumption of a linear relationship. The bars represent la error in 
fitting. 

The linear relationship given in Equation (1) is further tested by 
calculating the x 2 -, parametric Pearson's coefficient, rp, the rank cor- 
relation coefficient, rg, and their probabilities P p and P s respectively. 
Table II summarises the correlation statistics for all the data sets in- 
cluded in the fitting and shows that a positive correlation exists for 
all the activity indices. We further note that the best correlation is 
obtained for Fio, which confirms earlier results (Bachmann and Brown, 
1993; Bhatnagar, Jain, and Tripathy, 1999) that the radiative indices 
are better correlated with the frequency shifts. From the slope a and the 
gradient b obtained from the regression fitting, the following empirical 
relations between the shift, 5v and change in activity indices, Si are 
formulated: 
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Table II. Correlation statistics for Solar Cycle 22 
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5v = 


(2.44 ± 0.18) 5Ri - (6.00 ± 1.36) 


(2) 
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(2.41 ± 0.19) 5R S - (0.48 ± 1.68) 


(8) 



where 8v is given in nHz and the change in activity indices have their 
standard units. The second term in the parenthesis indicates la error. 
As mentioned previously, these relations are valid in the frequency 
range of 1500 to 3500 /iHz and £ between 5 to 99. 

These relations are used to estimate the mean frequency shifts for 
solar cycle 21 through 23 by taking the frequency data set GM2 as 
the reference point. As an example, the estimated frequency shift ob- 
tained from Equation (8) for the smoothed sunspot number is plotted 
in Figure 2 and yields a shift of 0.37 ± 0.03 fiHz between the minimum 
and maximum of the solar cycle 21. This value is consistent with the 
earlier results by Palle, Regulo, and Roca Cortes (1989) and Elsworth 
et al. (1990), who had obtained a shift of 0.44 ± 0.06 ^Hz and 0.46 ± 
0.06 //Hz respectively for the same period, using the Doppler velocity 
data of low degree {I < 3) modes. From Figure 2, we notice that a 
very good agreement exists between the calculated values of 5v and the 
observed GONG frequencies (diamonds), however a small difference for 
BBSO (triangles) and LOWL (squares) is noticed. This difference may 
be interpreted as due to the use of different spectral lines in BBSO and 
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Figure 2. Comparison of estimated and observed mean frequency shift for cycle 21 
through 23. The solid line is the mean shift as calculated from Equation (8). The 
observed shifts are represented by symbols; triangles for BBSO, squares for LOWL 
instrument and diamonds for GONG. The vertical line within the symbols represents 
the magnitude of la error in the frequncy shift. 

LOWL instruments (Howe, Komm, and Hill, 1999) or different data 
reduction techniques. 

3.1. Comparison between cycle 22 and 23 

To investigate the validity of empirical relations for different cycles, 
we calculated the fitting parameters a and b for cycle 23, using the 
available GONG data in the frequency range of 1500-3500 /iHz and I 
between 5 to 99. These parameters for solar cycle 22 and 23 are given 
in Table III and a reasonable agreeement is found for both the cycles. 

Therefore, we propose that the derived linear relations are indepen- 
dent of solar cycle and can be used to estimate the frequency shifts 
for past, present and future solar cycles. This is illustrated in Figure 3, 
wherein we have plotted the calculated frequency shifts using Equations 
(2) and (4) and the observed frequency shifts for cycle 23. The dashed 
and solid lines represent the estimated 5v, obtained from the relation 
for International sunspot number and 10.7 cm radio flux respectively, 
whereas the frequency shifts from the GONG data, with reference to 
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Figure 3. The estimated and observed frequency shifts for 1995-1998. The dashed 
line shows the shift, as estimated from the sunspot number (Equation (2)) and the 
solid line for the 10.7 cm (Equation (4)). The diamonds represent the observed 
frequency shifts from GONG data in the frequency range of 1500 and 3500 piHz 
and £ between 5 to 99, where the mean shift is computed by taking GM2 as the 
reference. 



Table III. Values of gradient a (in nHz/activity) and intercept b (in nHz) obtained from 
least square fitting for solar cycle 22 and 23 
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Figure 4- (a) Estimated 5u for cycle 23 from (a) the predicted smoothed sunspot 
number from SGD web site (solid line) (b) predicted smoothed sunspot number and 
(c) 10.7 cm radio flux from MSFC web site. The dashed lines show error in Sv due 
to the errors in predicted activity indices. 



GM2, are shown as triangles. It is clear that the observed frequency 
shifts for cycle 23 are in close agreement with those obtained from the 
derived relations for cycle 22. 

In Figure 4, we have plotted the estimated frequency shifts for 
the current solar cycle 23, using predicted smoothed sunspot num- 
bers and 10.7 cm flux. The solid line in Figure 4a shows the pre- 
dicted shift using R s , as listed in Solar Geophysical Data web page 
( btt p : / / www .ngdc.noaa/ gov /stp /st p . html ) . The dashed lines are the 
predicted errors due to error in K s . Each smoothed sunspot number 
represents the average of two adjacent 12-month running mean of 
monthly means. The predicted value of sunspot number is based on 
the actual value of smoothed sunspot number for February 1999, which 
utilises the values of the averaged monthly means from August 1998 
through August 1999. The predicted frequency shift for the solar cycle 
23 is estimated to be 265 ± 90 nHz corresponding to the predicted 
maximum sunspot number of 118.1 ± 35. In Figure 4b and 4c, we 
also show the estimated bv for the current solar cycle using the pre- 
dicted smoothed sunspot number and 10.7 cm radio flux obtained from 
http: / / wwwssl.msfc . nasa, gov /ssl/pad/ solar / predict . htm . 
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In summary, we have obtained empirical relations between the 
change in solar activity indices and the shift in p-mode frequencies 
for cycle 22. From these empirical relations, we have estimated p-mode 
frequency shift for the current cycle 23. It will be of interest to see 
how our predicted value agrees with the actual observed value during 
solar maximum in 2000. It is further shown that these relations are 
independent of the solar cycle and hence can be used to estimate the 
change in frequency for past, present and future epochs, if the solar 
activity index is known. 
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